mo = AOI::aoi_get("Mount Saint Helens") %>%
AOI::aoi_buffer(.5)
elev = elevatr::get_elev_raster(mo, z = 12)
writeRaster(elev, "/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif", overwrite = TRUE)
wbt_slope("/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif","/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_slope.tif")
wbt_aspect("/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif","/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_aspect.tif")
wbt_hillshade("/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif","/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_hillshade.tif")
elevt = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_slope.tif")
slop = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_slope.tif")
aspe = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_aspect.tif")
hisd = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_hillshade.tif")
plot(elevt, axe = F, box = F, col = viridis::viridis(256), main = "elevation")

plot(slop, axe = F, box = F, col = terrain.colors(256), main = "slope")

plot(aspe, axe = F, box = F, col = rainbow(8), main = "aspect")

plot(hisd, axe = F, box = F, col = gray.colors(256, alpha = .8), main = "hillshade")
